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We study the quadrature uncertainty of the quantum elliptical vortex (QEV) state using the 
associated Wigner function. Deviations from the minimum uncertainty states were observed due to 
the absence of the Gaussian nature. In our study of the entropy, we noticed that with increasing 
vorticity, entropy increases for both the modes. We further observed that, there exists an optimum 
value of ellipticity which gives rise to maximum entanglement of the two modes of the QEV states. 
A further increase in ellipticity reduces the entropy thereby resulting in a loss of information carrying 
capacity. We check the validity of the entropic inequality relations, namely the subaddivity and the 
Araki-Lieb inequality. The later was satisfied only for a very small range of the ellipticity of the 
vortex while the former seemed to be valid at all values. 
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I. INTRODUCTION 



Many of the fundamental concepts of quan- 
tum mechanics have been adapted from classical 
mechanics, albeit with some modifications. One 
such concept of particular importance is the phase 
space, which remains useful when passing to quan- 
tum mechanics. Similar to the probability den- 
sity distribution functions in classical mechanics, 
quasi-probability distributions were introduced in 
quantum mechanics [1 . Unlike the former, these 
can be negative which justifies the name, quasi- 
probability distribution functions. Among them, 
the Wigner function stands out because it is real, 
nonsingular and yields correct quantum mechani- 
cal operator averages in terms of phase-space inte- 
grals and possesses positive definite marginal dis- 
tributions [2_. The Wigner distribution function 
has come to play an ever increasing role in the 
description of both coherent and and partially co- 
herent beams and their passage through first order 
systems [3]. Once the Wigner function is known, 
other properties of the system can be calculated 
from it. 

Another concept, which quantum mechanics de- 
rives from classical mechanics, is that of entropy. 
It is a natural extension of the classical concept 
when dealing with quantum statistical mechanics. 
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It is also a key concept in the field of quantum in- 
formation theory. The entropy of quantum states, 
described in terms of density operators, replacing 
the classical probability theory, is defined by von 
Neumann entropy It gives a measure of how 
much uncertainty there is in the state of a physical 
system and also gives an idea about the informa- 
tion carrying capability of a state which is the task 
of present study. 

Optical vortices, possessing orbital angular mo- 
mentum (OAM), have been studied classically, by 
using phase plate or computer generated hologram, 
for quite some time [5] . Agarwal et al gave the con- 
cept of quantum optical vortex (QOV) , which were 
circular in shape [BHS]- It was later generalized to 
quantum elliptical vortex (QEV) [§]. In a recent 
work, it has been pointed out that photon sub- 
straction from one of the spontaneous parametric 
down converted beams (idler) produces an ellipti- 
cal vortex state [TU]. It should be highlighted that 
the photon subtraction from spontaneous paramet- 
ric down converted light has been realized experi- 
mentally as well |lll I12| . 

The paper is organized as follows. In section II 
we introduce the concept of QEV states in terms 
of basis vector states. We write the corresponding 
state for an elliptical vortex. We briefly discuss the 
Wigner distribution function for generalized quan- 
tum vortex. Using the same function, we calculate 
the uncertainty products and discuss the results 
therein. We write down the corresponding reduced 
density matrices for QEV states in section III and 
study the entropy of the constituent modes. In 
section IV we verify the validity of the entropic 
inequalities, namely, subadditivity and Araki-Lieb 
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inequality for the QEV states. We conclude this 
article after pointing out the significant results and 
directions for future work in section V. 



II. UNCERTAINTIES OF THE TWO 
MODES OF QUANTUM ELLIPTICAL 
VORTEX (QEV) STATES 

The Gaussian wave packets occupy a central 
place in studies involving wave packets of quantum 
system. In case of radiation fields these packets 
play an important role as these are in fact mini- 
mum uncertainty states and describe both, the co- 
herent states [TJ] as well as the squeezed states 
|14| . The quantum elliptical vortex states, as in 
the following equation, are of non-Gaussian struc- 
ture 




where m is an integer, measuring vorticity, or the 
maximum OAM of the state. 

The squeezed state \?p) for the two mode radiation 
field is defined as the direct product of the two 
squeezed mode states \i/j) a and \4')b, 



= exp [C,{a t2 - a 2 }] |0) 
x exp [C,{6 t2 - b 2 }] |0>. 



(2) 



Here a, b are the regular bosonic annihilation oper- 
ators for the two modes and Q's are the squeezing 
parameters. Invoking the disentangling theorem 
|14| . the exponential for one mode in Eq. Q can 
be expressed as 



exp [( x (a t2 - a 2 )] = exp 



^a t2 



exp j— ln[cosh(2£ x )] ( a^a + 



exp 



tanh.(2£ z 



(3) 



which is a product of exponentials. Using the fact 
that a and b acting on their respective vacuum 
states give zero, we can write the squeezed state 
for mode a as 



1 



= , . exp 
v /cosh(2C 2; ) 



t2 



10). (4) 



Following a similar approach we can also write 
the corresponding expression for \ip)b- The vortex 
state can now be written as 



Wqev) = A (jfeat - iriytf) 



(5) 



where A is the normalization constant, repre- 
sents the two squeezed mode state as in Eq. (§, 
m stands for the vorticity and rj x , rj y control the 
ellipticity of the vortex. If we put r\ x — rj y = 1, 
( x = Cy = C (real), it reduces to the circular 
vortex state (QOV) in a Gaussian beam. Using 
r/i = 1/ (v2crj), the normalized spatial distribu- 
tion of the QEV state is obtained as 



QEV 



2(m-2) 



x exp 



where <Ji = exp(2£. 




trarily and change the variables to a new set of 



scaled ones defined as P, 



P, 



and P„. 



^Py M- Following the 



treatments of [15 , the four dimensional Wigner 
function for the state ^>QEv{x,y), is obtained in a 
compact fashion as, 

W(x,y,p x ,py) 

= K exp [- (X 2 + Y 2 + P 2 1 + P 2 J] 

\2 



(P X2 + P Y2 -x 2 - Y 2 y 



(7) 



where K = ?T, [-2 (a\ + a 2 v )\ m and 



--1/2 



is the associated Laguerre polynomial. 



Using Eq. Q we determine the uncertainty in 




Figure 1: (Color online) Variation of uncertainties in 
quadratures in one mode with a x 

x, y, p x and p y . We study the uncertainties as a 
function of a x for different vorticities m. It is seen 
from Fig. [I] that A X increases monotonically with 
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Figure 2: (Color online) Variation of uncertainties in 
quadratures of the other mode with a x 



uncertainty relation. It is further observed that at 
a particular value of a x all the uncertainties are 
same. It is noticed from Fig. [3] that it is no longer 
the minimum uncertainty state. Hence the uncer- 
tainty product AxAp x decays to l/y/2. AyAp y 
has an almost complementary nature to that of 
the former. It starts from around l/y/2 and gradu- 
ally increases to reach saturation at a much higher 
value. The maximum value attained by AxAp x is 
different than that attained by AyAp y . The same 
holds true for the minimum values as well. In Fig. 
[2] we study the sum of the uncertainty products, 
AxAp x + AyApy. We observe that it nearly re- 
sembles a displaced Gaussian. 
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Figure 3: (Color online) Uncertainty products of both 
the modes with <j x 




Figure 4: (Color online) Sum of the uncertainty prod- 
ucts with a x 



increasing a x whereas Ap x starts from infinity and 
falls to a much lower value with increasing a x . This 
is expected as Ax varies linearly with a x and Ap x 
varies as l/(Ax) and hence l/(er x ). On the other 
hand Ay increases nonlinearly with a x (Fig. [2]). 
This is mainly because of the parametrization that 
we used for our calculations where Ay varies as 
y/5a x . Ap y on the other hand starts from infinity 
and falls as 1/yfdx, thus preserving the Heisenberg 



III. REDUCED STATE FOR EACH MODE 
AND INFORMATION ENTROPY 

In this section we calculate the reduced density 
matrices of the two modes and calculate the cor- 
responding entropies. We start with determining 
the expression for the density matrix of the QEV 
state \4>qev)- Using Eq. Q and (pj), one can write 
down the vortex state in an expanded form in Fock 
state basis, as 

\iPqev) = A(rj x a* -iriytf)™ '\i/j) a \ip)b 
A 

cosh(2Cx)cosh(2^ J/ 



E 

fe=0 



k\(m — fc)! 



{-my) Vz 



6c 

x ex P 1 y a 



t2 



|0)a 



ft , 



6 t2 ]6t fc |0> fc . 



x exp 

The density matrix p of the vortex state is 
P = \i/jqev)(^qev\- 



(8) 



(9) 



Using the properties of creation operators a) and 
and their conjugates we can write down Eq. ^ 
in the following form 

A 2 



cosh(2£ x ) cosh (2( y 

m ! 2 



X 



E " t! 2(m-fe) 2fc 



k=0 

x exp 
x exp 



fc!(m — fc)! 
^y« t2 j \m - k) a a (m- k\ exp 



tf* \k) 



b b\ 



(10) 



To calculate the entropy for each mode we need to 
obtain the corresponding reduced density matrices, 
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which is defined by p a = Tr^p). In this context we 
should mention that the reduced state is a mixed 
state even though the two mode state of Eq. Q 
is a pure state. Tracing out the b mode, we can 
write down p a as 



Pa 



^x 

cosh(2£i;)cosh(2(^ 



E 

fe=0 



k\(m — k)\ 



2(m-fe) 2fc 
'ix ly 



x exp \ —a 



,(k\ exp 



t2 ^j \m - k) a a (m - k\ exp ( 



£,x 2 



(ii) 



where A x is the normalization constant for mode a. 
Solving the term in square bracket and rearranging 
we obtain 



A 2 



Pa 



cosh(2Ca;)cosh(2C a ) 

m |2 

ml 



x 



2(m-fc) 2fc 

T77Z IXi^x >n 



k=o 
X exp 

x 2 Fi 



k\{m-k)\ lx 



l fc 



- i f 2 



(12) 



where 2 Fi [^r 1 , ^J 2 -, 1,£ 2 ] is the Hypergeometric 



function [IB]- Using Eq. (12 1 we calculate the 



diagonal elements of p a . We make use of the fact 
that, Trp a = /CfcLo w h ere are the diagonal 
elements of p a . 

On rearranging the terms we finally arrive at the 
following form 



Trp a 



E^ 

A;=0 
m 

E 



A 2 



k=0 



cosh(2( a; )cosh(2( a 



2(m-fc) 2fc 
x 77 77. 



» (m-fc)! 



E 

n=0 



^\ 2 " 1 



(to — k)\ 



n! 2 (to — fc — 2n)! 



x( to — fc — 2n\m — k) (to — k\m — k — 2n) 



2F1 



1 fc + 2 



1 f 2 



(13) 



Invoking the orthogonality criterion we see that 
only the n — term will remain while all other 



terms will cancel. So we can write down the sim- 



plified form of Eq. ( 13 ) 

m 

Trp Q = ]Tq? 

k=0 



A 2 



E cosh(2C*)cosh(2Cj,) 



2(m-fc) 2fc 
X 7 /x Vy 



TO ! 



X 2 Fl 



(to — /c)!fc! 

' 9 > x >Sy 



(14) 



Using these coefficients we calculate the von Neu- 
mann entropy of mode a, following the treatment 
of [7] as follows 



(15) 



fc=0 



where C? stands for the coefficients in ( 14 1. The 



logarithm is taken to base 2 as is the norm for in- 
formation entropy. 

Following a similar approach we can determine the 
von Neuman entropy for the b mode. Repeating 
the procedures of Eq. (Ill and Eq. (12) with the 



only difference that we trace out the a mode in- 
stead of &, we can write the reduced density matrix 
Pb = Tr Q (p) as 



Pb 



cosh ( 2 ) cosh ( 2 (j, ) 



E 



x exp 



x 2F1 



fc!(m — k) 



2(m-fe) 2k 
I 'lx 'ly 



tf 2 )\k) bb (k\ 



; l exp 



£1/1,2 



TO — fc + 1 TO— fc + 2 



l £ 2 



(16) 



where A y is the corresponding normalization con- 
stant for mode b. Using Eq. ( 16 ) we determine the 



diagonal elements of the reduced density matrix for 
the b mode as follows 



Trp 6 



E^ b 

m 

E 



A 2 



k=0 



cosh (2 £ x ) cosh (2£ y ) 

ml 2 



2(r 



-fc)„2ft 
'ly 



{m-k)\k\ 
m — k ■ 



x 2F1 



1 to - k + 2 



- 1 £ 2 



(17) 
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The corresponding entropy, Sb, can be calculated 
using Eq. ( 17 1 and the equivalent form of Eq. ( 15 I 



by replacing a with b. 
We study the corresponding entropies with re- 




Figure 5: (Color online) Variation of Entropy(S'a) with 



spect to rj x . We use a y — 3, a x = 5 in our cal- 
culations for generating the graphs. We use the 
same parameterizations for Q and £j as in section 
II, where the subscript i stands for x and y. We 
choose rj y — l/(^/2r] x ) arbitrarily. The normaliza- 
tion constant is evaluated for the two modes as 



A- 



{l + al)(l+a y ) 2 k\(m-k)\ 



(18) 



where the subscript i is used to denote the modes 
a and b. Fi represents the corresponding Hyperge- 
ometric function as expressed in Eq. ( 12 1 for mode 
a and in Eq. ( 16 1 for mode b. 



We have studied the entropies of the two modes 




Figure 6: (Color online) Variation of Entropy^) with 

as functions of r\ x in Fig. [5] and Fig. [6] for differ- 
ent m, the vorticity. We observed that the entropy 
increases with increasing vorticity. It is also ob- 
served that the peaks occurred at almost a fixed 
value for each mode for all values of m although it 



occurred for different values of r\ x for the modes a 
and b. This signifies that there exists an optimum 
value of rj x for which maximum entanglement can 
be achieved for the QEV state. In other words, 
an optimum level of ellipticity exists for which we 
can obtain the maximum entanglement. We use 
the word entanglement as information entropy is a 
measure of the degree of entanglement of the con- 
stituent states which make up the system of states. 
A further look at Fig. 5 and Fig. 6 makes us no- 
tice another very interesting fact. The maximum 
entropy occurs, not at the value of rj x = 1, but at 
some other value. It is to be noted that r) x = 1 
corresponds to the circular vortex, which has been 
a topic of study for long. This observation leads 
us to conclude that elliptical vortex has more en- 
tropy than the circular vortex which means more 
information transfer is possible by elliptical vor- 
tices, choosing correspondong ellipticity, than the 
circular vortex. This observation further empha- 
sises the importance of the need to study elliptical 
vortices. 

If the ellipticity is further increased, entropy falls 
off exponentially. Thus we can conclude that, in- 
creasing the elliptic nature of the vortex beyond 
a certain value would lead to disentanglement and 
hence loss of information carrying capacity. 



IV. ENTROPIC INEQUALITIES AND 
THEIR VALIDITY FOR QEV STATES 

In the previous section we have calculated the 
reduced density matrices for the two modes using 
the partial trace operation and also calculated the 
corresponding entropies. In this section we check 
the validity of the entropic uncertainty relations 
for the QEV state. 

If two systems A and B have a joint quantum state 
Pab , then the entropy of the combined states are 
expected to satisfy the following inequalities [4. 



S(A,B) < 
S(A,B) > 



S(A)- 
\S(A) 



■3(B) 
-S(B)\ 



(19) 
(20) 



We adapt these equations to suit our system where 
these are equally valid as the subsystems of the 
QEV state, p a and pi, are distinct quantum states 
though they are correlated. So we expect the Sub- 
additivity, Eq. (19 1 to hold with the inequality 
The second equation, i.e. Eq. (20 1 is the Araki- 
Lieb inequality. This inequality has some impor- 
tant implications for the index of correlation [18] . 
The index of correlation, I c , is a measure of the 
information content of the correlation between the 
components of an N component system. For a two 
component system it can be expressed in a simpli- 
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Figure 7: Variation of Entropy with r/ x for m = 1 



ficd form 



S a + Sb — Sab 



(21) 



where S^b stands for the entropy of the combined 
system. If the two component system is in a pure 
state, it is shown that maximum correlation oc- 
curs |191 EU], The Araki-Lieb inequality serves an 
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Figure 8: Variation of Entropy with rj x for m 



important purpose here by limiting the maximum 
value that I c can take (Eq. (21 1), thereby limiting 
the information content of the correlation of the 
components. 

To verify these inequalities we need to determine 
the entropy of the combined system that is, both 
the modes taken together. In section III we have 
written the density matrix of QEV state in Eq. 
(10 1. We use it to calculate the von Neumann en- 
tropy of the state. Proceeding as we did to deter- 



can find the trace of the entire system as follows 

m 

Trp = ]Tcf 



k=0 



A 2 



cosh(2( x )cosh(2('j / ) 
ml 



E 

fe=0 



(m — k)\k\ 



V 2 J m - k) V 2 y k (22) 



Eq. (22) is used to determine entropy of the 
combined system, i.e., S a b- We use the same 
parametrization of the previous section to study 
them in Figs. (7][9). It is observed that, S a + Sb 
starts from zero and increases to reach saturation 
or a plateau region and falls off to zero. The joint 
entropy S a b on the other hand increases with in- 
creasing r\ x until it reaches a maximum. Then it 
starts decreasing gradually to reach zero. \S a — Sb\ 
exhibits a strikingly different nature than the other 
two. It has a couple of peaks with the same max- 
imum value while the other two attain the maxi- 
mum value only once. \S a — Sb\ starts from zero 
and increases to attain the maximum. After which 
it falls off to zero before rising again to the same 
maximum value for a different value of r\ x . Then 
it gradually falls off to zero like S a + Sb and S a b- 
The occurence of the zero at the middle can be ex- 
plained very easily. It is clear from Fig. [5] and Fig. 
[6] that both S a and Sb attain the same value for 
that particular value of rj x due to which \S a — Sb\ 
goes to zero. 

A further look at Figs. (7|9l enables us to infer 
that the entropies satisfy the Subadditivity, Eq. 
(19), for all values of r\ x but not with an equal- 



ity as is expected for correlated systems. This sig- 
nifies that Eq. (19 1 holds for QEV states. This 
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mine the trace of the reduced density matrices, we 



Figure 9: Variation of Entropy with rj x for m = 5 

was also expected as the entropy for an entangled 
state should be less than that of the summation of 
the entropies of the constituent systems. On the 
other hand the Araki-Lieb inequality, Eq. (20 1, 
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was satisfied only in a very small range of values of 
r\ x . For all other values it is violated. We can ar- 
gue that since the vortex is bounded by a sharply 
peaked Gaussian distribution, it is in this region 
that Eq. p0| is violated as verified by Keitel and 
Wodkiewicz[21 . Since r\ x controls the ellipticity 
of the vortex it is observed that the Araki-Lieb 
inequality is valid only for a very short range of 
ellipticity. In other words there exists a limited 
range of ellipticity for which Araki-Lieb inequality 
holds. As the ellipticity increases S a b also increases 
and attains a maximum value where the inequal- 
ity holds but as the ellipticity increases further the 
combined entropy falls off and the inequality is vi- 
olated. 



V. CONCLUSION 

In this paper we have calculated the uncertainty 
products using the Wigner function of the QEV 
state. We noticed that the uncertainty product 
AxAp x attains a minimum value of -i. It has 
a maximum value of about 1.25. AyAp yi on the 
other hand, has an initial value of about 1 / It 
starts increasing gradually and saturates at nearly 
1.2. It can be argued that the presence of the vor- 
tex modifies the characteristics of this state and 
hence it is no more the minimum uncertainty state 
which it would have been if the vortex was not 
present. 

We have studied the von-Neumann entropy of the 
QEV states in terms of basis vector states. We 
found that the entropy was raised considerably 
with the increase in the vorticity of the states. It 



was noticed that the peaks for both the modes oc- 
cur at different values of r\ x where r\ x is a measure 
of the ellipticity of the vortex. But the peak value 
of the entropies for the two modes remained the 
same. It was further observed that there exists 
an optimum value of ellipticity which gives rise to 
maximum entanglement of the two modes of the 
QEV states. A further increase in ellipticity re- 
duces the entropy thereby resulting in a loss of in- 
formation carrying capacity. 

We checked and verified the entropic inequalities. 
We observed that the strong subadditivity was sat- 
isfied for all conditions. This is expected for any 
entangled state, as the combined entropy for such 
state would be always less than that of the con- 
stituent systems taken together. We noticed that 
Araki-Lieb inequality, an indicator of degree of en- 
tanglement, was violated in all the regions, except 
in a very narrow region of ellipticity values. 
Our results serve as a pointer for further investiga- 
tion and studies of quantum elliptical vortices as 
means of information transport. The study of de- 
coherence in this system needs to be pursued, for 
finding out the robustness of this correlated sys- 
tem. Investigation of quantum discord in this case 
will also be exciting, from both fundamental, as 
well as application point of view. 
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